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We study the geometrical entanglement of the XXZ chain in its critical regime. Recent numerical 
simulations [Q.-Q. Shi, R. Oriis, J. O. Fjaerestad and H.-Q Zhou, New J. Phys. 12, 025008 (2010)] 
indicate that it scales linearly with system size, and that the first subleading correction is constant, 
which was argued to be possibly universal. In this work, we confirm the universality of this num- 
ber, by relating it to the Affleck-Ludwig boundary entropy corresponding to a Neumann boundary 
condition for a free compactified field. We find that the subleading constant is a simple function 
of the conipactification radius, in agreement with the numerics. As a further check, we compute it 
exactly on the lattice at the XX point. We also discuss the case of the Ising chain in transverse 
field and show that the geometrical entanglement is related to the Affleck-Ludwig boundary entropy 
associated to a ferromagnetic boundary condition. 









o 



> 



o 
o 



X 



I. INTRODUCTION 

The geometrical entanglement (GE) is a measure of the 
multipartite entanglement in a wave-function, and quan- 
tifies pthe distance to the closest untangled (separable) 
state.Era Starting from a wave-function \^) for a quan- 
tum lattice model with N sites, one defines a maximal 
overlap A,„ax to be 



A, 



max|($|'i')|, 



(1) 



where the maximization is carried on states |$) which 
are tensor products of single-site states : 



N 



l<I'>=(g)|1'.) 



The larger A,„ax, the closest it is to a product state, p^nd 
the less entangled |^) is. The GE is then defined asE 



Em) = -\og.,k 



2 
max • 



This quantity allows to quantify the global multipartite 
entanglement of the wave function | ^) , and is useful in 
the context of quantum computationd or state discrimi- 
nation with local measurementsQ 

The GE has also recently gained, interest in many- 
body and condensed-matter physicsErQ. Since it gener- 
ally scales linearly with the volume [N] of the system 
one also defines the GE per site: 

£Nm) = N-^Em)- 

As an example of interesting application, it has been 
shown that derivatives of £m on large systems can be 
used to detect quantum phase transitionsJ3n3. 

In this paper we focus on the case of the ground-state 
JVP) = \G) of the Hamiltonian of a periodic spin-1/2 XXZ 
chain: 



H 



E(-J- 



j+i 



y V 



Aa|a|+i] 



(2) 



where L = N is the number of sites. We consider the 
anisotropy parameter in the range JA| < 1, so that the 
system is gapless (critical) . In Ref. ^ it was show numer- 
ically that the GE per spin admits the following asymp- 
totic expansion : 



Sl{A) = £^{A) + b{A)/L + 0{l/L^) 



(3) 



The purpose of this paper is to show that the term 6(A) 
appearing in this expansion can be computed analytically 
in a simple way. Upiig results of boundary Conformal 
Field Theory (CFT)lia, we relate it to an Affleck-Ludwig 
(AL) boundary entropy.E3 



II. RELATION TO THE AL BOUNDARY 
ENTROPY 

The product state |$)max that maximizes the overlap 
in Eq. |l| with the ground-state of Eq. || is kncpar to be 
a spin configuration where all spins are paralleltS and lie 
in the XY plane : 
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-L/2 



I free) J 
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Computing 6(A) amounts to calculating a subleading 
contribution to the scalar product ^(freejG). To do this, 
we adopt a transfer matrix point of view where \G) is 
interpreted as the dominant eigenstate of the transfer 
matrix M of some two-dimensional classical system (6- 
vertex like in our case). The classical model is defined 
on a cylinder of circumference L and height Ly ^ L, 
with free boundary conditions for the spins (or the 6- 
vertex arrows) at both edges. In this geometry the. 
free energy F = — ln^(frec|M'^y |free)z can be writtenEHl 
as F = i^buik + 2i^boundary, with Fbuik ~ LLy and 
^boundary = aL + s + o(l). o is the boundary free en- 
ergy per unit length, and s is a subleading term in the 
boundary free energy. 



It is easy to check that i^boundary = — In ^(frce|G) in the 
limit Ly ^ L where only the state \G) contributes. So, 
6(A) is simply related to the constant s in the boundary 
free energy: 



^(^)— h^^ 



(4) 



For critical systems s is universal and may therefore 
be computed in the continuum limit. In the XXZ chain 
(equivalently 6- vertex case) , it corresponds to a compact- 
ified free field with the following (Euclidian) Lagrangian 
density: 

an 
0=0 + 2ttR 

The compactification radius R is related to the decay ex- 
ponents of the correlation functions and to the Luttinger 
parameter. In the case of Eq. H, i? is a known function 
of AM 



i?(A) = W-arccos(A). 



(5) 



As well-known in boundary CFTIIj, free boundary con- 
dition for the spins (or arrows) are encoded in the con- 
tinuum limit by a Neumann boundary condition for the. 
free field. So, the boundary entropy s is in fact the ALt2l 
boundary entropy sn associated to Neumann boundary 
condition. Its value is knownlij and depends solely on the 
compactification radius: 



SN = ln(57v) , Qn = V^/2, 

where g^ is the so-called "universal non-integer ground- 
state degeneracy"J13 or "^-factor". Combining this with 
Eq. Q, our prediction for 6(A) is : 



6(A) = l-log2i?(A). 



(6) 



As can be seen in Fig. |i|, this result matches very well 
the numerical data of Ref. ^ The slight discrepancy for 
small A is likely due to finite-size effects and/or finite- 
bond dimension errors in the matrix-product represen- 
tation of the ground-state. This— can be confirmed by 
an exact microscopic calculationll2l of 6(A) at the free 
fermion point A = 0. 

The free fermkui representation of the chain at A = 
allows to obtainE£l £l on the lattice, as done in Ref. [Tp 
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From this one gets the expansion in powers of 1/L using 
the Euler-Maclaurin formula.Ej The result is: 



£l 



£. 

1- 
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6 

2K 

7rln2 



0(1/L2) 
~ 0.158733 



Subleading constant in the Geometrical entanglement for the XXZ chain 




FIG. 1. (color online) 6(A) for the XXZ chain, defined in 
Eq. H, as a function of the anisotropy parameter A. Symbols 
are numerical data of Ref. H while the solid curve is the field 
theoretical calculation (Eqs. m and pf). At A = 0, an exact 
lattice calculation shows that 6 = 1. 



where K is Catalan's constant. The result can be com- 
pared with the numerical estimations of Ref. M : E^o — 
0.1593 and 6 ~ 0.98. As expected, the exact lattice cal- 
culation confirms the field theory prediction at A = 0, 
namely 6(A = 0) = 1 (Eq. |). 

It is also worth noticing that the calculation extends to 
finite magnetizations of the chain, and simply amounts 
to move the fermion density away from 1/2. In such a 
case one observes that 6(A = 0) = 1 is independent of 
the magnetization. This again matches the field theory 
result, since, for A = 0, the compactification radius R is 
indeed independent of the magnetization. 



III. CRITICAL ISING CHAIN 

The GE for the periodic Ising chain in transverse field 
at the critical point has also been considered numerically 
in Ref. Bl The Hamiltonian is now 



^ = -E-.^ 



L 



(8) 



and the product state which majtimizes the overlap turns 
out to be a tilted configuration :I3 



N 



l*)r 



i(cos(^/2)|t,)+sin(^/2)|;,)) (9) 



C=i 0.897101 



where | fj) and | ]^j) are the eigenstates of cr^.Ej The 
authors of Ref. ^jfound numerically in this case 6 ~ 1.016. 
We will now give an argument based on boundary CFT, 
similar to that presented in the case of the XXZ chain, 
which shows that 6=1. 



Again, we wish to interpret the scalar product A,„ax 
appearing in the GE as a boundary contribution to a 
classical two-dimensional free energy. When \G) is the 
ground state of Eq. ^, the corresponding classical model 
is a critical two-dimensional Ising model. If we had to 
project \G) onto a state where all spins would point in 
the X direction (corresponding to an angle ^ = ±7r) , it 
would correspond to a fixed ferromagnetic boundary con- 
dition for the Ising model. On the other hand, if we had 
to project onto a state with all the spins pointing in the z 
direction (^ = 0), it would correspond to a free boundary 
condition (as in the case of the XXZ chain) . We therefore 
see that the tilted state of Eq. somewhat corresponds 
to a combination of free and fixed boundary conditions 
for the classical model. But the important point is that 
^ 7^ 0. For this reason. |<E>)max breaks the Z2 symmetry 
of the model, which exchanges the x and —a; directions. 
In such situation, where the boundary condition imposes 
a non-zero magnetization at the edge, the long-distance 
and universal properties of the boundary will be equiva- 
lent to that of a system with fixed ferromagnetic bcuind- 
ary condition (all spins pointing in the x direction) J13 In 
other words, the constant term s in the boundary free 
energy will be the same for the tilted boundary condi- 
tion and for the ferromagnetic one. In this case, the AL 
entropy Sfixcd is known to be - 
).00 
we therefore obtain b = 1. This argument can also be 



' 2 In 2 (corresponding to 
a g-factor gnxed = l/V2)-t3'Eil As conjectured in Ref. ^, 



confirmed using the exact formula of Ref. M for ±he GE, 
and then applying the Euler-Maclaurin formula.E3 



dimensional Rokhsar-Kivelson wave functions, as studied 
in Refs. ^ and gi] (the Renyi parameter being n = (3/2). 
Using numerical and field-theoretical approaches, it was 
found that In Z{(3) scales as In Z{(3) = aL + ^, where the 
constant term 7 is universal and given by 



7=(l-/3/2)lni?+-ln(/3/2). 



(11) 



So, for the XXZ chain, the GE appears to be a special 
case /3 = 1 of this partition function. 

We finally comment on the the limit /3 — )■ 00. where 
only the spin configuration |i)max of the z basis which 
have the largest weight in \G) contributes to Z{/3). For 
the present XXZ chain, |i)max is the ferro mag netic state 
|i)max = I t • • • t> -1^3 Taking /3 ^ 00 in Eq. 0and Eq. |ll| 
we get ln(t • • • t \G) = — ^Ini? (omitting the extensive 
term). So, we see that the scalar product of \G) with both 
(i) the uniform configuration where spins point in the z 
direction and (ii) the configuration where spins point in 
the X direction, have very similar subleading terms re- 
lated to the logarithm of the compactification radius. In 
the GET language, the former scalar product (appearing 
in Z{l3 = 00)) corresponds to a Dirichlet boundary con- 
dition for the free field, as opposed to Neumann for the 
scalar product appearing in the GE and 2'(/3 = 1). 



CONCLUSION 



IV. RELATION WITH RENYI ENTROPIES FOR 
ROKHSAR-KIVELSON IN 2+1 DIMENSIONS 

As we have seen in Sec. ||, the GE for a critical XXZ 
chain can be expressed using a sum of scalar products 
with all the spin configurations \i) of the z-basis: 



Z = 2^/2 A 



E(*|G)- 



In conclusion, we have related subleading terms of the 
geometrical entanglement of the critical XXZ and Ising 
spin chains to their AL boundary entropies, showing that 
this correction is universal. Whether such universal be- 
haviour in the scaling of the geometrical entanglement 
(or of the scalar product of the ground-state with a given 
separable state in general) can be found for other physical 
systems is an open issue that is worth pursuing. 



Z can be seen as a partition function and generalized by 
introducing a fictitious inverse temperature /? > 0£3 



Z(/3)=E((*|G))^ 



(10) 



It turns out that Z{f3) is related to Renyi entangle- 
ment entropies of a semi-infinite cylinder in some two- 
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